Less singular quasicrystals: life in low codimensions 
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We consider a set of tilings proposed recently as d- 
dimensional generalisations of the Fibonacci chain, by Vidal 
and Mosseri. These tilings have a particularly simple the- 
oretical description, making them appealing candidates for 
analytical solutions for electronic properties. Given their self- 
similar geometry, one could expect that the tightbinding spec- 
tra of these tilings might possess the characteristically singu- 
lar features of well known quasiperiodic systems such as the 
Penrose or the octagonal tilings. We show here, by a numeri- 
cal study of statistical properties of the tight-binding spectra 
that these tilings fall rather in an intermediate category be- 
tween the crystal and the quasicrystal, i.e. in a class of almost 
integrable models. This is certainly a consequence of the low 
codimension of the new tilings. 

INTRODUCTION 

We report the results of a numerical study of the two- 
dimensional (2D) members of a new family of quasiperi- 
odic tilings proposed recently by Vidal and Mosseri |Q. 
These new so-called generalized Rauzy tilings (GRT) (af- 
ter the work of 0), are generalizations of the Fibonacci 
chain and are simpler in structure than hitherto stud- 
ied tilings in two and three dimensions. Tight-binding 
models on these systems may be amenable to analyti- 
cal techniques whereas other 2D tilings thus far studied 
have been refractory to approaches including renormal- 
ization group, or recursive methods ft0]. This property 
makes the new tilings attractive candidates to study in 
the world of quasiperiodic Hamiltonians. It is thus in- 
teresting to examine whether the tight-binding spectrum 
of these new tilings have similarly complex structures as 
for the "canonical" cases studied earlier, namely the Pen- 
rose or the Ammann-Beenkker (octagonal) tilings. This 
turns out not to be the case: these tilings fall in a class 
intermediate between the periodic lattice, an integrable 
problem, and the octagonal tiling, which belongs in the 
non-integrable class. 

This conclusion is based on calculations of statistical 
properties of the energy levels of periodic approximants of 
these new systems. Such analyses have been carried out 
for complex systems, including disordered media, systems 
with interactions or quantum billiards Q. A relation is 
found to exist between the quantum statistics and the 
classical dynamics : depending on whether the dynamics 
is classically ergodic, or classically chaotic, one finds their 
quantum properties to be described by Poisson statistics 
or by Random Matrix Theory (RMT) . Such statistical 



analyses are useful for gaining insight into complex sys- 
tems for which exact, or even approximate theoretical re- 
sults are not available. This is the case for the octagonal 
tiling, a quasiperiodic tiling in two dimensions which has 
been tackled via numerical study of tight-binding mod- 
els 0-^]. We show that the GRT in two dimensions has 
spectral properties that differ essentially from those of 
the octagonal tiling. The GRT appears analogous to the 
pseudointegrable billiards proposed in |Tc| which are nei- 
ther completely integrable nor completely chaotic. 

DESCRIPTION OF THE NEW TILINGS 

The cut-and-project method provides a simple and 
easily implemented way to obtain pieces of quasiperi- 
odic tilings. One can thus obtain arbitrarily large pieces 
of an infinite quasicrystal, or one can construct finite 
pieces that can be periodically continued - these are 
the periodic approximants of the quasicrystal. As de- 
scribed in ipH , one first defines the real-space directions, 
which constitute a d-dimensional subspace of a larger D- 
dimensional hyperspace. The remaining D-d directions 
constitute the 'perpendicular" space. The quasicrystal is 
obtained by projecting all the points contained within a 
certain infinite strip in the hyperspace Jl^j onto the real 
space. 

When the orientation of the real space is irrational, the 
tiling is aperiodic, while if rational orientations are cho- 
sen, the tiling repeats periodically in real space. This 
method has been used to obtain the Fibonacci chain 
(D=2, d=l), the octagonal tiling (D=4, d=2), the three- 
dimensional Penrose tiling (D=6, d=3) etc. The codi- 
mension, D-d, is in fact equal to d in the three cited 
cases. 

The simplifying feature of the GRT is that their 
codimension is equal to 1 for all real space dimension 
d. For a finite tiling, this allows one to index the points 
according to their perpendicular space coordinate. Fur- 
thermore, when the points of an approximant of these 
tilings - illustrated for two dimensions in Fig.l - are so 
indexed, the connectivity matrix of the approximant (as- 
suming toroidal boundary conditions) takes on a very 
simple band-diagonal (Toeplitz) structure 0. In the two 
dimensional case that we consider here, the nonzero ma- 
trix elements are situated at distances of F n -3, F-n—2 and 
F„_i from the diagonal. The generalized Fibonacci num- 
bers {F n } are obtained from a three term recursion rela- 
tion F n = F„_i + F„_2 + -Fn-37 with the initial con- 
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ditions F_i = 0; F = Fi = 1. The ratio F„/F„_i 
tends to the value a w 1.839..., solution of the equa- 
tion x 3 — x 2 + x + 1. These matrices correspond to 
approximants of increasing size as n is increased. 

The model we consider is a pure hopping tight-binding 
Hamiltonian, 

H = Yi( it j)tcJ Cj + h.c. (1) 

where i,j in the sum correspond to pairs of sites linked 
by a bond in Fig.l. The hopping amplitude has been as- 
sumed to be independent of the linked sites and is set to 
unity. We consider the case of periodic boundary condi- 
tions, and take k = where k is the Bloch vector. The 
resulting eigenvalue problem reduces to that of diagonal- 
izing the connectivity matrix defined above. This was 
done using a Lanczos routine. The spectrum is expected 
to be symmetric for the infinite tiling, since it is bipar- 
tite. In addition, there is a discrete symmetry under an 
inversion in the perpendicular subspace |l3| . This results 
in two subspaces of positive and negative parity respec- 
tively, and our analyses are carried out for each sector 
separately. 

RESULTS 

The spectral density p(E) is shown in Fig. 2a, with 
that of the octagonal tiling shown in 2b for compari- 
son (the localized states of the latter at E = have 
not been included in the figure). The density of states 
curves shown were obtained for given sample sizes assum- 
ing periodic boundary conditions, and illustrate the char- 
acteristic rapidly fluctuating behavior typical of these 
aperiodic tilings (similar curves were found in Q| and 
0). It is evident that while the GRT has an underlying 
smoothly varying component of p{E) which is everywhere 
non-zero (notice the y-axis offset) whereas the octagonal 
tiling does not. Calculations on approximants of the 2D 
Penrose tiling show strong fluctuations as well with a 
multitude of gaps and pseudogaps [|l4|. 

The envelope of p{E) of the GRT is in fact reminiscent 
of the density of states of the square lattice. One notes 
in Fig.l that there are many regions in which sites of 
coordination number 4 are grouped together in this tiling. 
In contrast to the octagonal tiling, which also has a mean 
coordination number ~z = 4, the fluctuations of geometry 
are less strong in the GRT. This is a result of the reduced 
codimcnsion, and we see that in consequence the density 
of states has a much less singular structure as compared 
to the octagonal case. 

Turning to the statistics of nearest neighbor level spac- 
ings, we investigate the distribution of E i+ i — Ei where 
the Ei are the ordered set of energy levels of the Hamil- 
tonian (1). In order to eliminate the trivial dependence 
of this quantity on variations of p(E), one first carries 



out an "unfolding" of the spectrum ]l5[ to get the cor- 
rected level spacings [|6| . This yields the unfolded spac- 
ings Si whose mean value is equal to 1. The resulting 
distribution of spacings P(s) is plotted in Fig. 3a for the 
three largest systems studied, which contain 35738 (L), 
19513(M) and 10609(S) sites. The dashed line repre- 
sents the Poisson e _s decay, while the continuous line is 
the semi- Poisson form Psp{s) = 4sexp(— 2s). The latter 
has been proposed for pseudointegrable billiards and can 
be shown to correspond to a certain short range plasma 
model (SRPM) by Bogomolny et al JTH]. These authors 
speculate that this type of statistics is a consequence 
of fractal structure of the wave functions in those sys- 
tems. In the Anderson model, the critical level statistics 
were first shown to have the semi-Poisson distribution by 
Braun et al |]l7f . In another example of this universality 
class Evangelou and Pichard had considered the distribu- 
tion of bandwidths and obtained the semi-Poisson form 
in the critical Harper model fl8[| . 

In Fig. 3a one sees that there is a size-dependence in 
the fall-off in P(s) at small s. For the largest size, the 
points appear to deviate substantially from the semi- 
Poisson curve. It will be interesting to calculate the next 
system size to see if the deviations are systematically 
upwards, i.e. towards the Poisson limit. The P(s) re- 
sembles that obtained in Jl7| for periodic boundary con- 
ditions at the critical point of the Anderson model. In 
contrast to the random case however, in our case there 
is a size-dependence. It is important to note that the 
degree of unfolding matters here: nonunfolded spacings 
follow a simple Poisson law while progressively stronger 
unfolding causes the distribution to go over to the semi- 
Poisson form. In Fig. 3b we present a semilog plot of 
the statistics of the maximally unfolded level spacings ( 
where we renormalized each spacing by the locally aver- 
aged value over the nearest neighbor energy levels) . The 
agreement with the semi-Poisson form is seen to be good 
over a wide range of values. At large s P(s) falls off 
faster than exponentially in s, while at small s, there is 
a size dependent increase of P(s) above the semi-Poisson 
value. We believe that these effects may be explained 
by the unfolding procedure, although boundary condi- 
tion effects cannot be ruled out and further studies are 
needed to settle this issue. 

This distribution in Fig. 3a is very different from the 
one obtained for the octagonal tiling. For that tiling, 
unfolded spacings follow the Wigner-Dyson distribution 
for random matrices of the gaussian orthogonal ensemble 
(GOE), Pgoe(s) = | exp(-7rs 2 /4). If one considers the 
distribution of the "bare" spacings (i.e. without correct- 
ing for the variation in density of states), they follow a 
broad (log-normal) distribution law : this can be shown 
via a recursive relation based on the inflation symme- 
try of the tilings p9| . A less singular case is that of a 
randomized version of the octagonal tiling Is. The ran- 
domized tiling has a less fluctuating density of states and 
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the spacings distribution is simply GOE (unfolding does 
not produce any new effects). From these P(s), as well 
as from the spectral rigidity function (to be discussed be- 
low) we conclude that the octagonal tilings belong in the 
class of non-integrable models exemplified by disordered 
metals. One has an unbounded diffusion of wave packets 
resulting in level repulsion at small s. The GRT is clearly 
not in the GOE regime for any scale of energies, rather it 
should be considered a nearly integrable system, as the 
shape of P(s) indicates. 

This becomes more evident upon looking at the spec- 
tral rigidity function £( 2 )(L) = ((N(L) -N) 2 }, where 
N(L) is the number of levels lying within an energy in- 
terval of width L. The angular brackets mean averages 
over different starting positions of the interval within the 
spectrum. Fig. 4a shows the plot of T,^ for the largest 
tiling size, for three different choices of unfolding |T(| 
The continuous line in Fig. 4a indicates the SRPM result 
£( 2 )(L) = i(L + (1 - e~ 4L )/4). For very small L, the 
dependence appears well described by the semi-Poisson 
law. The three curves have the same low energy behav- 
ior, but the "most unfolded" case has its fluctuations 
damped out because of the unfolding, while the other 
cases show the transition from approximately linear to 
approximately quadratic dependence on L. 

The SRPM behavior at small L may be related to a 
long time dynamics in this tiling which is intermediate 
between chaotic and integrable. This would correspond 
to the semi-Poisson form of P(s) ft is interesting to con- 
sider the rigidity evaluated for the energy levels without 
unfolding. In the case of random media, where diagram- 
matic perturbation theory can be exploited, and checked 
by more semiclassical methods |2(| , it was shown that in 
the diffusive regime, as a function of energy E, T,^(E) 
has two different sorts of behavior: firstly, at very low 
energies, it has the logarithmic dependence expected for 
GOE matrices. This holds upto the scale of the Thouless 
energy E c pf|] which is size dependent, corresponding to 
the time taken to diffuse to the boundaries. For larger 
energies S^ 2 * 1 ~ E du where the exponent v describes the 
spreading over time of an initially localized wave packet 
(r 2 ) ~ t 2v . This theoretical prediction was verified nu- 
merically in the Anderson model [p2[ , where the exponent 
is one for ordinary diffusion, v — 0.5. For the random- 
ized octagonal tiling, v w 0.85 J^] as deduced from E^ 2 \ 
in agreement with the results of a direct calculation of 
wavepacket diffusion on this tiling [^3). Fig. 4b shows 
the log-log plot of YP^ as a function of energy for three 
sizes. There is a crossover from linear dependence to an 
almost quadratic dependence, YP^{E) tx E 2v {v rs 0.95). 
The sub-quadratic regime very likely corresponds to the 
sub-ballistic regime of dynamics that was observed for 
wavepackets diffusing on a 2D GRT by Vidal and Mosseri 




DISCUSSION AND CONCLUSIONS 

The two-dimensional GRT tiling, despite its quasiperi- 
odic structure, is an almost integrable system, with a 
semi-Poisson form for the level spacings distribution and 
a corresponding small-energy spectral rigidity function. 
The P(s) distribution is probably linked to the existence 
of fractal wavefunctions, as has already been observed in 
critical quantum systems. In the GRT studied here, it 
is noteworthy that although fluctuations of the density 
of states are present, they are strongly reduced with re- 
spect to those in "canonical" quasiperiodic tilings. This 
is a consequence of the reduced codimension of the tilings, 
which leads to its having weaker quasiperiodic modula- 
tions of geometry than the previously studied octagonal 
and Penrose tilings. Further study is needed to elucidate 
the nature of states and of the relation between dynami- 
cal evolution of wavepackets and the calculated statistical 
properties of the energy levels. 
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FIG. 1. An approximant of the 2D rauzy tiling showing the 
numbering scheme: the arrows show the diplacement vector 
from site 0( origin) to site 1, etc. Note the periodic boundary 
conditions. 
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FIG. 2. (a) Density of states of the 2D Rauzy tiling and 
(b) of the octagonal tiling 
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FIG. 3. (a) P(s) of three approximants: L(circles), 
M (diamonds) and S (triangles) along with the semi-Poisson 
(continuous line) and the Poisson (dashed line) laws, (b) 
Semilog plots of P(S) with the Poisson and the semi-Poisson 
laws (dashed and continuous lines respectively) . 
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FIG. 4. (a) The rigidity E 2 (£) for different levels of unfold- 
ing, (see text) from most unfolded (circles) to not unfolded 
(triangles), (b) E 2 (_E) without unfolding for three sizes L,M 
and S. The dashed lines show E and E 2 power laws. 
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